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Transport of fermions is central in many fields of physics. Electron transport runs modern tech- 
nology, defining states of matter such as superconductors and insulators, and electron spin, rather 
than charge, is being explored as a new carrier of information J^. Neutrino transport energizes 
supernova explosions following the collapse of a dying star [2], and hydrodynamic transport of the 
quark-gluon plasma governed the expansion of the early Universe However, our understanding 
of non-equilibrium dynamics in such strongly interacting fermionic matter is still limited. Ultracold 
gases of fermionic atoms realize a pristine model for such systems and can be studied in real time 
with the precision of atomic physics [4] [5]. It has been established that even above the superfiuid 
transition such gases fiow as an almost perfect fluid with very low viscosity [31 [5] when interactions 
are tuned to a scattering resonance. However, here we show that spin currents, as opposed to 
mass currents, are maximally damped, and that interactions can be strong enough to reverse spin 
currents, with opposite spin components reflecting off each other. We determine the spin drag coeffi- 
cient, the spin diffusivity, and the spin susceptibility, as a function of temperature on resonance and 
show that they obey universal laws at high temperatures. At low temperatures, the spin diffusivity 
approaches a minimum value set by h/m, the quantum limit of diffusion, where h is the reduced 
Planck's constant and m the atomic mass. For repulsive interactions, our measurements appear to 
exclude a metastable ferromagnetic state [THS] . 



Understanding the transport of spin, as opposed to 
the transport of charge, is of high interest for the novel 
field of spintronics [1] . While charge currents are unaf- 
fected by electron-electron scattering due to momentum 
conservation, spin currents will intrinsically damp due to 
collisions between opposite spin electrons, as their rela- 
tive momentum is not conserved. This phenomenon is 
known as spin drag [lOl [H] . It is expected to contribute 
significantly to the damping of spin currents in doped 
semiconductors [T^. The random collision events also 
lead to spin diffusion, the tendency for spin currents to 
fiow such as to even out spatial gradients in the spin den- 
sity, which has been studied in high-temperature super- 
conductors US] and in liquid ^He-^He solutions [TS] . 

Creating spin currents poses a major challenge in elec- 
tronic systems where mobile spins are scattered by their 
environment and by each other. However, in ultracold 
atoms we have the freedom to first prepare an essen- 
tially non-interacting spin mixture, separate atoms spa- 
tially via magnetic field gradients, and only then induce 
strong interactions. Past observations of spin currents in 
ultracold Fermi gases jTH [17] were made in the weakly- 
interacting regime. Here we access the regime near a 
Feshbach resonance [5], where interactions are as strong 
as allowed by quantum mechanics (the unitarity limit). 
We measure spin transport properties, the spin drag coef- 
ficient Fsd and the spin diffusivity Dg, of a strongly inter- 
acting Fermi gas composed of an equal number of atoms 
in two different spin states. In the strongly-interacting 
regime, spin drag is expected to reach a universal max- 
imum value, and spin diffusion is expected to reach a 
universal minimum. 

The universal behaviour of the spin transport coeffi- 



cients of a Fermi gas can be estimated on general grounds. 
At the Feshbach resonance, the scattering cross-section 
a between atoms of opposite spin is given by the square 
of the matter wavelength, in the degenerate regime a ^ 
l/fc|,, where kp = (Ott^tt,)^/'^ is the Fermi wavevector and 
n is the density of atoms in each spin state. The mean 
free path between collisions is thus / = 1/na ^ l/kp 
or about one interparticle spacing, the smallest possible 
mean free path in a gas. The average speed v of atoms is 
on the order of the Fermi velocity Kkp /rn. In estimating 
the spin diffusivity Dg « vl the density-dependent factors 
cancel, giving Ds ~ h/m. This value for represents 
a universal quantum limit to spin diffusivity in Fermi 
gases. Away from resonance, the scattering cross sec- 
tion decreases, increasing Dg. For temperatures T much 
greater than the Fermi temperature Tp — h^kp/2mkB, 
the scattering cross section will be given by the square 
of the thermal de Broglie wavelength and thus decrease 
as a cx 1/r, while the velocity will increase as w oc 
causing Dg to increase as Dg cx T3/2. Finally, in a degen- 
erate Fermi gas the average velocity will remain on the 
order of the Fermi velocity, but the effective scattering 
cross section will scale as ct oc due to Pauli blocking, 
causing to increase like as the temperature is low- 
ered. For a Fermi gas, we thus expect the minimum Dg to 
occur at temperatures near the Fermi temperature. Cor- 
respondingly, the coefficient Fsd characterizing spin drag 
is expected to reach a universal maximum value on reso- 
nance and for temperatures near the Fermi temperature, 
given by the Fermi rate Ep/h, where Ep — h?kp/2m is 
the Fermi energy. 

In the experiment, we prepare an equal mixture of 
the two lowest hyperfine states ( "spin up" and "spin 
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FIG. 1. Observation spin current reversal in a resonant 
collision between two oppositely spin-polarized clouds 
of fermions. (a) shows the total column density and (b) the 
difference in column densities of the two clouds (red: spin up, 
blue: spin down), in 1 ms intervals during the first 20 ms 
after the magnetic field is set to the Feshbach resonance at 
834 G. The collision leads to the formation of a high-density 
interface between the two spin states, (c) The separation 
between the centers of mass of the two spin states initially 
oscillates at a frequency of 1.63(2) oj^, where cJz = 27rx22.8 Hz 
is the trap frequency in the axial direction (see Sup. Info, for 
discussion). Even after half a second, there is still substantial 
spin separation. The diffusion time indicates a diffusivity on 
the order of h/m. (d) Shows the harmonic trapping potential 
along the axis of symmetry. 



down" ) of fermionic ^Li in a cylindrically symmetric atom 
trap [31 [5] . The confinement along the axis of symme- 
try is harmonic, with frequency lOz- We separate the 
two spin components along the axis of symmetry of the 
trap (see Methods Summary) and turn on strong inter- 
actions between unequal spins by quickly ramping the 
magnetic field to a Feshbach resonance located at 834 G. 
The confining potential of the trap forces the two clouds 
of opposite-spin atoms to propagate towards each other. 



establishing a spin current. Measurements are made by 
selectively imaging the two spin components. 

Figure 1 shows the collision between the two spin do- 
mains on resonance. The clouds bounce off each other 
and essentially completely repel. Due to the axial trap- 
ping potential, the clouds return after the collision, and 
we observe several oscillations in the displacement d — 
(z^) — (zi), where (^^(^)) is the center of mass of the 
spin up (down) cloud. After the oscillations have de- 
cayed, the displacement decreases to zero monotonically, 
on a timescale on the order of one second, an extremely 
long time compared to the trapping period (44 ms). The 
underlying explanation for spin current reversal and the 
slow relaxation can be found in the extremely short mean 
free path and the high collision rate between opposite- 
spin atoms at unitarity. According to the above estimate, 
the spin diffusivity is approximately h/m, which for ^Li 
is (100 /im)^/s. The atom clouds in the experiment have 
a length on the order of 100 /im, and it takes them on 
the order of a second to diffuse through each other. So 
we are indeed observing quantum-limited spin diffusion. 
The initial bounces will occur when the mean free path 
of a spin up atom in the spin down cloud is smaller than 
the spin down cloud size, i.e. when the mixture is hydro- 
dynamic. Instead of quickly diffusing into the spin down 
region, it is then more likely that the spin up atom is 
scattered back into the spin up region where it can prop- 
agate ballistically. Indeed, we see bounces occur already 
for interaction strengths away from the Feshbach reso- 
nance where the mean free path is on the order of the 
axial cloud size (see Supplementary Information). 

The relaxation dynamics close to equilibrium give di- 
rect access to the spin transport coefficients. The spin 
drag coefficient Tsd is defined as the rate of momentum 
transfer between opposite spin atoms |10j . and is there- 
fore related to the collision rate. The relaxation of the 
displacement d near equilibrium then follows the differ- 
ential equation [18^ 

Tsdd + uj'id = 0. 

Fitting an exponential with decay time r to the displace- 
ment gives the spin drag coefficient as Fgd — ^It. 

The spin drag coefficient is found to be greatest on 
resonance, and thus spin conduction is slowest (see Sup- 
plementary Information). On dimensional grounds, Tgd 
must be given by a function of the reduced tempera- 
ture T/Tp times the Fermi rate Ep/h. At high tem- 
peratures, we expect the spin drag coefficient to obey 
a universal scaling Tgd ~ nav ^ T^^/^. In fig. 2 we 
show the spin drag coefficient as a function of T/Tp. 
We observe T'^/^ scafing for T/Tp > 2, finding F^d = 
0.16(l)^(r/rF)"^/^- At lower temperatures, we ob- 
serve a crossover from classical to non-classical behav- 
ior: the spin drag coefficient reaches a maximum of ap- 
proximately O.lEp/h near the Fermi temperature. We 
interpret this saturation of the spin drag coefficient as 
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a consequence of Fermi statistics, as a and v approach 
constant values [16] determined by the Fermi wavevector 
kp- The maximum spin drag coefficient corresponds to 
a minimum spin conductivity Gg = on the order 

of kp/h. This is the slowest spin conduction possible in 
three dimensions in the absence of localization. 
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FIG. 2. Spin drag coefficient of a trapped Fermi gas. 

Spin drag coefficient for resonant interactions, as a function of 
the dimensionless temperature. The temperature is normal- 
ized by the maximum Fermi temperature, defined in terms of 
the maximum density of a single spin state in the system. We 
find agreement between measurements taken at three different 
values of the axial trapping frequency ujz — 2ir x 22.8 Hz (red 
circles), lj^ = 2tt x 37.5 Hz (blue triangles), ljz = 2^ x ff.2Hz 
(black squares). The data for T/Tp > 2 fits to a T"^/^ law 
(solid line). Dashed line: a power law fit with an uncon- 
strained exponent for T/Tp < 0.5 to show the trend. The 
spin drag coefficient is obtained from the decay rate of the 
spin separation and normalized by the square of the trap fre- 
quency and by the maximum Fermi energy. The error bars 
are statistical. 

At low temperatures, the spin drag coefficient de- 
creases with decreasing temperature. No theoretical pre- 
diction for the spin drag coefficient exists in this regime. 
For highly-imbalanced spin populations, a decrease in 
the spin drag coefficient is expected due to Pauli block- 
ing ^19J, which reduces the collision rate in the gas by a 
factor of [T/TpY at sufficiently low temperatures. How- 
ever, for systems with equal spin populations, it was 
shown that pairing correlations enhance the effective col- 
lision rate for collective excitations, leading to a collision 
rate that increases dramatically as the temperature is 
lowered [20] . The observed reduction of T^^ at low tem- 
peratures contrasts with this prediction, and is qualita- 
tively consistent with the onset of Pauli blocking. 

Comparing the relaxation rate to the gradient in spin 
density allows us to also measure the spin diffusivity D^- 
At the center of the trap, where the trapping forces van- 
ish, the spin current density Js is given by the spin dif- 
fusion equation [21] 



where n^{\.) is the density of spin up(down) atoms. 
We calculate the trap-averaged spin current as 
Js — \{n'^ + %)c?, where the densities are evaluated at 
the trap center. 

We find that spin diffusivity is minimum when interac- 
tions are resonant (see Supplementary Information). The 
increase in spin diffusivity for positive scattering length 
a, as well as the decrease in spin drag, argues against 
the existence of a ferromagnetic state in repulsive Fermi 
gases, for which diffusion should stop entirely [71 1^ I22j. 
Figure 3 reports the measured spin diffusivity as a func- 
tion of temperature at unitarity. In the high-temperature 
limit on resonance, one expects Dg ~ v/na oc T^/"^. At 
high temperatures, we indeed find this temperature de- 
pendence, with a fit giving = b.S{2)^{T /Tpf/"^ for 
T /Tp > 2. Down to our lowest temperatures, the spin 
diffusivity is seen to attain a limiting value of 6.3(3)/i/to. 
Fermi liquid theory predicts a rapid increase at low tem- 
peratures as Ds oc {Tp/TY'. From measurements on spin 
quadrupolar oscillations in highly spin-imbalanced mix- 
tures it is known [53j that the resulting crossover from 
diffusive to ballistic motion occurs below T/Tp k, 0.1. 
However, below T/Tp w 0.17 a balanced gas is superfluid 
and spin transport possibly inhibited by the superfluid to 
normal interface |21j . 
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FIG. 3. Spin diffusivity of a trapped Fermi gas. Spin dif- 
fusivity on resonance (solid circles) as a function of the dimen- 
sionless temperature T/Tp. At high temperatures, the spin 
diffusivity obeys the universal T-^/^ behaviour (solid line). At 
low temperatures the diffusivity approaches a constant value 
of 6.3(3)fi/m for temperatures below about 0.5rF, establish- 
ing the quantum limit of spin diffusion for strongly interacting 
Fermi gases. 

When comparing these results to theoretical calcula- 
tions, it is important to account for the inhomogeneous 
density distributions and velocity profiles that result from 
the trapping potential. For a homogeneous system on 
resonance, and at high temperatures compared to the 
Fermi temperature, we predict = 1.11^(T/Tf)^/^ 
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and r^d = 0.90^ {T/Tf)-^^^ (see Supplementary Infor- 
mation) . The measured spin drag coefficient is smaller by 
a factor of 0.90/0.16(1) = 5.6(4) while the spin diffusivity 
is larger by about the same factor, 5.8(2)/l.ll = 5.3(2). 
This factor reflects the inhomogeneity of the system and 
agrees with an estimate from the Boltzmann transport 
equation (see Supplementary Information). 

Finally, the measured transport coefficients give for the 
first time access to the temperature-dependence of the 
spin susceptibility Xs(^) in strongly interacting Fermi 



gases. Defined as Xs 



^("t-"4.) 



the spin susceptibility 



9(/i^-/i4,)' 

describes the spin response to an infinitesimal effective 
magnetic field or chemical potential difference fi^ — /i^ 
applied to the gas, and is a crucial quantity that can 
discriminate between different states of matter [5]. In 
a magnetic field gradient, particles with opposite spin 
are forced apart at a rate determined by the spin con- 
ductivity cTg, while diffusion acts to recombine them. 
The balance between the processes of diffusion and con- 
duction therefore determines the resulting magnetization 
gradient, a connection expressed in the Einstein relation 
Xs = o's/Ds [9]. In calculating this ratio from observable 
quantities, the relaxation time r cancels, as both as and 
are proportional to - , yielding 



Xs 



1 din-f - ni) 



'.di 



dz 



where 



d{r: 



is evaluated near the trap center. The 



inhomogeneous trapping potential has practically no ef- 
fect on the measurement of Xs, as all quantities involved 
refer to the vicinity of the center of mass. 

Figure 4 reports our findings for the spin susceptibil- 
ity at unitarity, as a function of the dimensionless tem- 
perature T/Tp. At high temperatures, we observe the 
Curie law Xs = n/T. In this classical regime of uncor- 
related spins, the susceptibility equals the (normalized) 
compressibility of the gas k = dn/dfj, that we also di- 
rectly obtain from our profiles. 

At degenerate temperatures, the measured suscepti- 
bility becomes smaller than the compressibility. This is 
expected for a Fermi liquid, where Xs — i+f" 
K = i^ps iSj with Landau parameters Fq and Fq 
describing the density (s) and spin (a) response. The 
spin susceptibility is expected to strongly decrease at suf- 
ficiently low temperatures in the superfluid regime, as 
pairs will form that will not break in the presence of an 
infinitesimal magnetic field. For the temperature range 
of our experiment this behaviour is not yet observable. 
It is currently under debate whether the strongly inter- 
acting Fermi gas above the superfluid transition temper- 
ature is a Fermi liquid |24) or a state with an excitation 
gap (pseudo-gap) [25l [26] . The opening of a gap in the 
excitation spectrum would be revealed as a downturn of 
the spin susceptibility below a certain temperature. Such 
a downturn is not observed in Xs down to T /Tp sa 0.2, 
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FIG. 4. Spin susceptibility on resonance determined 
from the Einstein relation, (a) Compressibility (solid blue 
circles) and spin susceptibility (open red circles) normalized 
by the compressibility of an ideal Fermi gas at zero tem- 
perature. For temperatures below the Fermi temperature, the 
susceptibihty becomes suppressed relative to the compress- 
ibihty, due to interactions between opposite-spin atoms. The 
spin susceptibility coincidentally matches the compressibility 
of a non-interacting Fermi gas (dashed line) in the range of 
temperatures that we could access, (b) Red circles: spin sus- 
ceptibility divided by the compressibility obtained from the 
same clouds. At temperatures above the Fermi temperature, 
the ratio of spin susceptibihty to compressibility approaches 
unity (dashed line). 



and therefore our spin susceptibility data agree down to 
this point with the expected behavior for a Fermi liquid. 

In conclusion, we have studied spin transport in 
strongly interacting Fermi gases. The spin diffusivity 
was found to attain a limiting value of about 6.3/i/m, 
establishing the quantum limit of diffusion for strongly 
interacting Fermi gases. Away from resonance the 
diffusivity increases. This casts doubt on the possibility 
of stabilizing a ferromagnetic gas on the repulsive side 
of the Feshbach resonance [7], which would require a 
vanishing diffusivity 0| . The observed slow relaxation of 
spin excitations is a likely explanation for the surprising 
- possibly non-equilibrium |21| - profiles in imbalanced 
Fermi gases observed at Rice [27], which did not agree 
with equilibrium measurements at MIT |28l [29] and at 
the ENS 23J. Our measurements of the temperature 
dependence of the spin diffusivity at low temperatures 
disagree with the expected behavior of a Fermi liquid. 
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while the spin susceptibihty that we measure is consis- 
tent with a Fermi liquid picture. An interesting subject 
of further study is whether spins are still able to diffuse 
through the superfluid, or whether they travel around 
it, avoiding the superfluid due to the pairing gap. 

Methods Summary 

The spin mixture is initially prepared at 300 G. To 
separate the spin components, we ramp the total mag- 
netic field to 50 G, where the magnetic moments of the 
two spin states are unequal, and apply two magnetic field 
gradient pulses. We then bring the total magnetic field 
to the Feshbach resonance in about 2 ms. 

To reach low temperatures during the approach to 
equilibrium, evaporative cooling is applied, at 834 G, by 
gradually lowering the depth of the optical dipole trap. 
To reach high temperatures, we heat the atoms by switch- 
ing off the optical dipole trap for up to 3 ms to allow the 
atoms to expand before recapturing them. We then set 
the final depth of the dipole trap so that the atom num- 
ber and the temperature remain nearly constant during 
the approach to equilibrium. 

Spin selective imaging is performed via in situ absorp- 
tion or phase contrast imaging using two 4 fj,s imaging 
pulses separated by 6 ^s. These images give the col- 
umn densities of each spin state, from which we obtain 
the three-dimensional density via an inverse Abel trans- 
form [29]. The gradient in the spin density is obtained 
from a linear fit to the polarization versus z. 

We determine the temperature of the clouds by fitting 
the density versus potential energy in the vicinity of 
z — 0, but for all r, to the equation of state of the unitary 
Fermi gas, measured recently by our group [3D]. The 
trapping potential itself is determined by summing the 
densities of hundreds of clouds, using the known axial, 
harmonic trapping potential to convert equidensity lines 
to equipotential lines and fitting the result to an analytic 
model. 
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DEPENDENCE ON INTERACTION STRENGTH 



We study the dependence of the spin transport properties of the system on interaction strength by ramping to a 
variable final field in the vicinity of the Feshbach resonance and measuring the subsequent evolution of the system. 
Figure 5 shows the results of colliding two clouds at different fields, revealing the transition from transmission of the 
clouds through each other to reflection of the clouds as the mean free path becomes smaller than the cloud size. When 
the scattering length is set to zero (Fig. 5. a), the center of mass separation oscillates at the trap frequency uj^. On 
resonance (Fig. 5.g), the observed oscillation frequency of 1.63(2) is intermediate between the frequency 1.55 uj^ of 
the axial breathing mode of a unitary Fermi gas in the hydrodynamic limit [5] and the non-interacting value of 2 , 
as the system contains a hydrodynamic region at the center, and is non-interacting in the spin-polarized wings. Figure 
6 shows the spin transport coefficients versus interaction strength. The spin drag coefficient exhibits a maximum on 
resonance (6. a), while the spin diffusion coefficient is minimum on resonance (Fig. 6.b). 
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FIG. 5. Collision between spin up and spin down clouds with varying interaction strength. After separating the 
spin components, the magnetic field was ramped to a variable value near the Feshbach resonance to reach different interaction 
strengths. The interaction parameter kpa, with kp = (Gn'^n)^^^ and n the central density per spin component, was determined 
by averaging the values of hp obtained from images taken after 200 ms of evolution time (not shown). The values of kpa were 
(a) 0, (b) 0.08, (c) 0.13, (d) 0.19, (e) 0.26, (f) 1.2, (g) oo, and (h) -1.5 




FIG. 6. Spin drag coefficient (a) and spin diffusivity (b) across the Feshbach resonance. The spin-separated clouds were cooled 
at 834 G to T/Tp ~ 0.16 (blue triangles) and to T/Tp ~ 0.40 (red cirlces) before ramping to the final field. The largest spin 
drag and smallest spin diffusivity occur at the Feshbach resonance, where 1/kpa — 0. The solid line in (a) is a Lorenztian fit. 
To reduce clutter, the 0.4 Tp data are not shown in (a). 
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IMAGE ANALYSIS 



The experiment takes place in a trapping potential of the form 

V{r, z) ~ -muj^z^ + Vr{r), 

where r = \/ +2/^- Here z is the symmetry axis of the trap, and we image along the y axis. 

From each run of the experiment we obtain two-dimensional column densities rt^^ix^ z) of both spin states a =t,4.- 
Fitting a two-dimensional Gaussian to each column density provides a measurement of the center of mass of each spin 
state. We subtract the z components of the centers of mass of the two spin states to obtain the separation parameter 
d = (z^) — {z^}. The three-dimensional densities n^(r, z) are obtained using a numerical implementation of the inverse 
Abel transformation (similar to [3T|): 



na{ri, z) 



2d 



ixj+i,z) - ni'^{xj,z) 



In 



j + - 



For each line of constant z, the values of n^iji^z^ with small i are sensitive to noise in n'^ {xj,z) for small j. To 
reduce noise, we fit a one-dimensional Gaussian to n'^'^{xj, z), for each value of z and for j < 4 (corresponding to the 
innermost 5 /im of the clouds or less than a tenth of the Fermi-radius), and use this fit (with sub- resolution sampling 
to reduce discretization error) in the above formula. 

The temperature is determined by analyzing the densities n„{r,z) for \z\ < 60 fxni (as the 1/e radii of the clouds 
range from 100 to 400 /im, this restriction excludes the more polarized, large z, regions of the clouds). We obtain the 
temperature of each cloud by fitting the density versus potential energy V to 



where (3 = jp^, A 



27r 



^j^^rp , T is the temperature, n is the chemical potential, and / is a universal function defining 
the equation of state of the unitary Fermi gas [30], with fj, and T as the fit parameters. 



MEASUREMENT OF TRANSPORT COEFFICIENTS 



Near equilibrium, d decays exponentially, and we fit the measured values to 

d{t) = doe"*/^. 

The spin drag coefficient is then 

To non-dimensionalize this quantity we multiply by h and divide by the Fermi energy Ep = ■^{STT'^n{0))'^^'^ , where 
n{0) — n^{0) + "^(0) is the total three-dimensional density at (r, z) = (0, 0). 

After fitting to d{t) to obtain r, we obtain the spin diffusivity from each spin up-spin down image pair as 



n{6)d 
2gT ' 



where g — ^^"g^"''"'^ \z=o,r=o is the spin density gradient. The spin density gradient is obtained from the slope 6 of a 
linear fit to the polarization p{z) = ^^^^^ as a function of z at r = using the formula 

5 = 5-(n^(0)+n;(0)). 
The Einstein relation provides the spin susceptibility as 

n(0) 



Xs = (7s/ Ds = 



2m Ds Fsd 



idu!^ 



Although derived using transport coefficients, the final expression for Xs does not actually depend on the relaxation 
time r, so it gives an independent value of Xs from each run of the experiment. 
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THEORETICAL CALCULATION OF TRANSPORT COEFFICIENTS 



The transport coefficients in tire classical limit T 3> Tp may be calculated using the Boltzmann transport equation. 
The Boltzmann transport equation describes the evolution of the semi-classical distribution functions w, t), with 

— f + V ■ Vrfa H • V^/o- = /cou[/t4 ] 

ot m 

Define the joint distribution function /(f^, f;, w^, iT;) = f^{r-t,v^)fi{ri,vi). Any quantity of the form x{r^,ri,v^,vi) 
can be averaged over / [32j: 

(X) = ^ / d'r^d\id%d\(xf) , 

where is the number of atoms with spin a. 

We generalize to arbitrary scattering length a, and non-uniform drift velocity, the calculation by Vichi and 
Stringari [18^ of the equation of motion for the relative coordinates of the two spin components in a harmonic trap of 
angular frequency Uz along the z axis. From the Boltzmann equation, 

and 

dt {vzt - Vzl) = - {vzt - Vzl) Fsd - {z^ - zi) , 

where the spin drag coefficient is 

^ {{Vzt - Vzi)Icoll[f]) 

J-sd = ] r ■ 

The full collision integral reads 

[(1 - /t)(i - - /t/i(i - /;>(! -/!)]■ 

The integral describes the scattering of particles t and J, into new states f' and I'. We assume that the distribution 
functions near equilibrium are fair, v) — f'^{r, v — Va-.di-m), where /" is the equilibrium distribution of both spin states 
(assumed to be in the classical limit) and Wcr, drift = ivdi'r)z is the drift velocity of spin a. The collision integral can 
then be calculated (for details, see for example [33]) and yields 

Tsd = ^a(r/r,) (^) ' no(0)/a, 
6 \ irm J 



with 

a 

no = J /"d'^w is the equilibrium density of each spin state, and Ta = ^^^^^2 is the temperature scale associated with 
the scattering length a. The limiting values of the temperature-dependent average scattering cross section a{T/Ta) 
are 



= Atto^ du "'^ 



1 

1 -I- 4^u^ 



J noVd d^r 
j rinWd d'^r 
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and 

mKBi 

with the thermal de Broghe-wavelength A. For unitary interactions, Ta — 0, and the latter limit applies. Define the 
Fermi energy as Ep = ^(67r^no(0))^/'^. In the unitary limit where cr = we find 

hV^A _ 32^/2 _ 0.90 Fk^ 

1^ ~ 97r3/2aV A^bT ^ ^yksf' 

For a uniform system a = 1, while for a harmonically trapped system with a uniform drift velocity we find a = 2?!'^ . 
However, the drift velocity profile cannot be uniform: Even if it started out uniform, spin currents would get damped 
faster in the center of the overlap region of the two clouds, where the collision rate is high, than in the wings, 
where it is low. A non-uniform drift velocity profile will develop. The nature of the linearized Boltzmann equation 
allows for a variational principle where trial functions replacing the true distribution function j„ (r, v) yield upper 
bounds on the actual spin drag coefficient [331 . Minimizing Fgd for a trial class of non-uniform drift velocities 
va(v) — Wc;o(l — no(r)'''/no(0)''') yields 7 — 0. In this case a = 2^/^ 5.7. 
The spin diffusivity measured in our experiment for T ^ Tp is then 



1 no(0) 



32^2 m\Tp 
= 1.1a 



h f 3/2 



TO \Tp 

The bulk value for a = 1 represents a lower bound on the diffusivity. It is known from similar calculations that 
in the high-temperature limit, this lower bound should be within a few percent of the actual bulk value 133]. The 
effect of inhomogeneities in a harmonic, cylindrically symmetric trap increases the diffusivity by a w 5.7. The 
high-temperature result for the diffusivity has been obtained independently by Georg Bruun 34 . 

I 



